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ABSTRACT 

Here we report on the results of deep X-ray follow-up observations of four 
unidentified 7-ray sources detected by the Fermi /LAT instrument at high Galac- 
tic latitudes using the X-ray Imaging Spectrometers on-board the Suzaku satel- 
lite. All of the studied objects were detected with high significance during the 
first 3-months of Fermi /LAT operation, and subsequently better localized in the 
first Fermi/LAT catalog (1FGL). For some of them, possible associations with 
pulsars and active galaxies have subsequently been discussed, and our observa- 
tions provide an important contribution to this debate. In particular, a bright 
X-ray point source has been found within the 95% confidence error circle of 
1FGL J1231. 1-1410. The X-ray spectrum of the discovered Suzaku counterpart 
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of 1FGL J1231. 1-1410 is well fitted by a blackbody with an additional power-law 
component. This supports the recently claimed identification of this source with 
a millisecond pulsar PSRJ123 1-1411. For the remaining three Fermi objects, 
on the other hand, the performed X-ray observations are less conclusive. In the 
case of 1FGL J1311. 7-3429, two bright X-ray point sources were found within the 
LAT 95% error circle. Even though the X-ray spectral and variability properties 
for these sources were robustly assessed, their physical nature and relationship 
with the 7-ray source remain uncertain. Similarly, we found several weak X-ray 
sources in the field of 1FGL J1333. 2+5056, one coinciding with the high-redshift 
blazar CLASS J1333+5057. We argue that the available data are consistent with 
the physical association between these two objects, although the large positional 
uncertainty of the 7-ray source hinders a robust identification. Finally, we have 
detected an X-ray point source in the vicinity of 1FGL J2017.3+0603. This Fermi 
object was recently suggested to be associated with a newly discovered millisec- 
ond radio pulsar PSR J2017+0603, because of the spatial coincidence and the 
detection of the 7-ray pulsations in the light curve of 1FGL J2017.3+0603. In- 
terestingly, we have detected the X-ray counterpart of the high-redshift blazar 
CLASS J2017+0603, located within the error circle of the 7-ray source, while we 
were only able to determine an X-ray flux upper limit at the pulsar position. All 
in all, our studies indicate that while a significant fraction of unidentified high 
Galactic latitude 7-ray sources is related to the pulsar and blazar phenomena, 
associations with other classes of astrophysical objects are still valid options. 

9 Subject headings: galaxies: active — pulsars: general — radiation mechanisms: non- 
10 thermal — gamma-rays: general — X-rays: general 



11 1. Introduction 

12 Observations with the EGRET instrument onboard the Compton Gamma-Ray Observatory 

13 (CGRO) in the 1990's opened a new window in studying MeV-GeV emissions from both Galactic 

14 and extragalactic objects. Despite over a decade of multi- wavelength fol low-up studies, more than 



15 60% of the 7-ray emitters included in the 3rd EGRET catalog (3EG; iHartman et al.l I1999I) are 

16 yet to be identified (that is, 170 out of 271). This is mainly because of the relatively poor 7-ray 

17 localizations of EGRET sources (typical 95% confidence radii, r 95 ~ 0.4° — 0.7°), challenging the 
is identification procedure especially for the objects located within the Galactic plane, due to source 

19 confusion. In particular, as much as ~ 90% of the 3EG sources detected at Galactic latitudes 

20 |6| < 10° do not have robustly selected counterparts at lower frequencies. On the other hand, most 
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21 of the 3EG sources at high Galactic latitudes have been associated with blazars — a sub-class of 

22 jetted active galactic nuclei (AGN) displaying strong relativistic beaming — in accordance with 



23 the exp ectation for the extragalactic population to dominate the 7-ray sky at |6| > 10° (|Abdo et al 



2009 a|). Yet the unidentified fractio n of the high Galactic latitude 3EG sources is still large (~ 30%; 



25 e.g., 



Sowards-Em merd et al. 200 3J). The s ituation is basically unchanged in the revised EGRET 
26 catalog (EGR: ICasandjian & Grenier 2008 ). even though the revised background modeling applied 



27 in the EGR resulted in fewer 7-ray detections (188 sources in total, in contrast to 27 1 listed in 3EG); 

28 87 out of 188 EGR entries remain unidentified. 

29 The unidentified low Galactic latitude 7-ray sources are expected to be associated with lo- 

30 cal objects such as molecular clouds, supernova remnants, massi ve stars, pulsars and pulsar wind 

31 nebulae, or X-ray binaries (see, e.g., iGehrels & Michelson|[l999l and references therein). Mean- 

32 while, the population of unidentified high Galactic latitude 7-ray sources is typically believed to 

33 be pr edominantly extragalactic in origin, although there is a suspected Galactic component as 

34 well (lOzel & Thompsonlll996|) . For example, the brightest steady source 3EG J1 835+5918 lo- 



35 cated at \b\ > 10° was proposed to be associated with an isolated neutron star (iMirabal et al. 



2000j; iReimer et al.ll200l|) . The neutron star origin and its association with the 7-ray source has 



37 been confir med with the disco very of a 7-ray pulsar at the position of 3EG J 1835+59 18 with 

38 Fermi/L AT (A bdo et al .1 120 1 Oallbh . Sim ilarly, high-energy 7-ray pulsations were discovered with 

39 Fermi (kbdo et alJl2009bh and AGILE dTavani et al.ll2009h fro m PSR J202 1+3651 that was long 

40 considered as a likely pulsar counterpart of 3EG J2021+3716 (|Halpern et al.ll2008l) . On the other 

41 hand, blazar G74. 87+1.22 (B 2013+370) was claimed to be the most lik ely counterpart of the 

42 unidentified object 3EG J2016+3657 located within the Galactic plane (|Mukherjee et al.l 12000 : 



Halpern et al. 



2001) . Other unidentified 7 -rays sources were similarly investigated with varying 



Muk herjee & HalpernlEooH) . We note that population studies, which could in princi 



44 success (e.g., 

45 pie shed some light on the galactic/extragalactic origin of different classes of unidentified EGRET 

46 sources, were impeded by different level of background emission at different locations from the 

47 Galactic plan e, and different EGRET exposure for various parts of the sky (see the discussion in 



48 



Reimenl200l|) . Also, variability stu dies were previous ly hampered by the limited statistics and 



49 noncontinuous EGRET observations (INolan et al.l l2003). 



50 



With the successful launch of the Fermi Gamma-ray Space Telescope, we now have a new op- 

51 portunity to study 7-ray emission from different types of high energy sources with much improved 

52 sensitivity and localization capabilities than with EGRET. With its field of view (five-times-larger 

53 than that of EGRET) covering 20% of the sky at every moment, and its improved sensitivity 

54 (by more than an o rder of magnitude with respect to EGRET), the Large Area Telescope (LAT; 



ss lAtwood et al.l l2009) aboard Fermi surveys the entire sky each day down to a photon flux lev- 
seels of -F>iooMeV — few xl0~ 7 phcm~ 2 s _1 . The first FermifLAT point source catalog (1FGL) 
57 already surpasses EGRET with 1451 sources detected at significance levels > 4a within the 
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58 100 MeV — 100 GeV photon energy range during the initial 1 1 -month survey (|Abdo et al.ll2010ch . 

59 Several high-latitude EGRET sources lacking low-frequency counte rparts were confirm ed by Fermi/LAT 

60 and associated with previously unknown 7-ray blazars, as expected ( Abdo et aDboiOd ). Somewhat 

61 surprisingly, however, a number of 7-ray emitters at | b\ > 10° have be en robustly identifie d by LAT 

62 with newly found 7-ray pulsars via the detection of 7-ray pulsations (lAbdo et al.ll2010el) . Most of 

63 these are in fact millisecond pulsars (MSPs). A diminishing, yet still significant population of 

64 unidentified Fermi/LAT objects remains, constituting as much as about 40% of all 1FGL sources. 

65 This includes more than 10 unidentified EGRET sources at high Galactic latitudes, which are thus 

66 the best candidates for the persistent, or even "steady" 7-ray emitters over the 10-year-long period 

67 between the EGRET and FermifLAT epochs (as indicated by their comparable photon fluxes in the 

68 3EG and 1FGL catalogs). 

69 Thus motivated, we started a new project to investigate the nature of unidentified high Galactic 

70 latitude Fermi objects through deep X-ray follow-up observations with the Japanese X-ray astron- 



71 omy satellite Suzaku (|Mitsuda et al.ll2007|) . This paper presents the results of the first year cam- 

72 paign conducted over the span of Suzaku-AOA (Apr 2009 - Mar 2010), during which we have ob- 

73 served four s teady/weakly varia ble FermifLAT sources from the 3-month Fermi/LAT Bright Source 

74 List (0FGL; lAbdo et al]|2009c|) . These are denoted below accordingly to their 1FGL catalog en- 

75 tries as 1FGLJ1231. 1-1410, 1FGLJ131 1.7-3429, 1FGL J 1333. 2+5056, and 1FGLJ20 17. 3+0603. 

76 Thanks to the superb localization provided by the LAT, all the corresponding 95% error cir- 
77cles (typically r 95 ~ 0.1° — 0.2°) could be covered within the field-of-view of the Suzaku X- 

78 ray CCD camera "XIS". Only in the case of 1FGL J1333. 2+5056, the Suzaku pointing does 

79 not cover the entire 95% LAT error circle since the localization error for this object did not 
so improve sufficiently between 1FGL and 0FGL. Along with our Suzaku observations, system- 
si atic pulsar searches with radio telescopes have been performed for the Fermi/LAT unassociated 
82 sources. These resulted in the new discoveries of MSPs co-located with the two 7-ray sources 
ss included in our study (1FGLJ1231. 1-1410 and 1FGL J2017. 3+0603). In both cases, Fermi/LAT 
84 eventually detected 7-ray pulsations as w ell, in accordance with the results in the radio domain 
ss (|Ransom et al.ll2010l : ICognard et al.ll2010|) . Our deep X-ray exposure discussed in the next sec- 

86 tions supports the pulsar identification for at least 1FGL 1231.1-1410, but is less conclusive in the 

87 case of 1FGLJ2017. 3+0603. For the other target from our list, 1FGLJ1333. 2+5056, a tentative 
ss association with b lazar CLASS J 1333+5057 was claimed in the LAT Bright AGN Sample (LB AS; 
8 9 lAbdo et alJb009ah . Here we substantiate this possibility by presenting the broad-band spectral en- 



90 ergy distribution (SED) for 1FGL J1333.2+5056/CLASS J1333+5057, including new Suzaku data, 

91 which is indeed typical of a flat spectrum radio quasar (FSRQ). Finally, the nature of the remain- 

92 ing source 1FGLJ131 1.7-3429 (for which no radio or 7-ray pulsations have been detected so far; 

93 Ransom et al. 2010) could not be revealed, despite the discovery of a likely X-ray counterpart. In 

94 particular, we found that the multiwavelength spectrum of 1FGL J13 11.7-3429 is not consistent 
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95 with neither a typical blazar nor pulsar spectrum. 

96 In § 2, we describe the Suzaku X-ray follow-up observations and the data reduction procedure. 

97 The results of the analysis are given in § 3. The discussion and conclusions are presented in § 4 

98 and §5, respectively. A standard ACDM cosmology with = 0.73, Vl M = 0.27, and H = 

99 71 km s -1 Mpc -1 is assumed throughout the paper. 



100 2. Observations and Analysis 

101 2.1. Observations and Data Reduction 



102 
103 
104 
105 
106 
107 
108 
109 



We observed four unidentified high G alactic latitude Fermi/LAtI objects with the Suzaku X- 
ray astronomy satellite (|Mitsuda et al. Il2007h . These are denoted in the 1FGL catalo g as 1FGLJ123 1 . 1- 
1410, 1FGLJ131 1.7-3429, 1FGLJ1333. 2+5056, and 1FGL J2017. 3+0603 (see kbdo et ap2010ch . 
All the sources but one (1FGLJ20 17. 3+0603) were already listed in the 3rd EGRET catalog 
(IHartman et al.l 1 1999Q and their 7-ray fluxes are given in Table [Q The Suzaku observation logs 
are summarized in Tabled The obs ervations were made with three out of four CCD cameras 
(X-ray Imaging Spectrometers; XIS: iKovama et al.l |2007|) . and a Hard X-ray Detector (HXD; 



Kokubun et a l. 2007; iTakahashi et al.ll2007|) . One of the XIS sensors is a back-illuminated CCD 



1 
1 
1 
1 
1 

120 
121 
122 

123 
124 



(BI; XIS1), and the other three XIS sensors are front-illuminated ones (FI; XIS0, XIS2, and XIS3; 
the operation of XIS2 has been terminated in November 2006). Since none of the studied sources 
have been detected with the HXD, in this paper we focus on the analysis of only the XIS data. 
The XIS was operated in the pointing source mode and in the normal clocking mode during all the 
exposures. 

In the reduction and the analysis of the Suzaku data, HEADAS software version 6.7 and a 
calibration database (CALDB; released on 2009 September 25th) were used. The XIS cleaned 
event dataset was obtained in the combined 3x3 and 5x5 edit modes using xselect. We 
excluded the data collected during the time and up to 60 seconds after Suzaku was passing the South 
Atlantic Anomaly (SAA). We also excluded the data corresponding to less than 5 degrees of the 
angle between the Earth's limb and the pointing direction (the Elevation Angle; ELV). Moreover, 
we excluded time windows during which the spacecraft was passing through the low Cut-Off 
Rigidity (COR) of below 6 GV. Finally, we removed hot and flickering pixels (using sisclean; 



Day et al.lll998l) . With all the aforementioned data selection criteria applied, the resulting total 



effective exposures for all the observed sources are summarized in Tabled 



-6- 



125 



2.2. Analysis 



126 XIS images for each target were extracted from the two FI CCDs (XISO, XIS3) within the 

127 photon energy range from 0.4 to lOkeV. In the image analysis procedure, calibration sources lo- 

128 cated at the corners of CCD chips were excluded. The image s of Non X-ray B ackground (NXB) 

129 were obtained from the night Earth data using xisnxbgen (|Tawa et al.ll2008l) . Since the expo- 

130 sure times for the original data were different from that of NXB, we calculated the appr opriate 

131 exposure-corrected original and NXB maps using xisexpmapgen (|Ishisaki et alj 120071) . The 

132 corrected NXB images were next subtra cted from the correc ted original images. In addition, we 

133 simulated flat sky images using xi s s im (|Ishisaki et al. I l2007h . and applied a vignetting correction. 

134 All the images obtained with XISO and XIS3 were combined and re-binned by a factor of 4. All 

135 the FI XIS images were in addition smoothed by a Gaussian function with o = 0/17, and the 

136 resultant images are presented in section [3l Note that the apparent features at the edge of these 

137 exposure corrected images are undoubtedly spurious due to low exposure in those regions. For the 

138 further analysis, source regions were carefully selected around each detected X-ray sources within 

139 the error circle of a studied 7-ray emitter. The corresponding background regions with radius 3' 

140 were taken from the same XIS chips avoiding any bright X-ray spots in the field. In all the cases, 
hi such source regions were set to within 3' or V radii around the X-ray point sources (because of 

142 the blurring due to the Suzaku/XIS Point Spread Function; PSF), depending on the properties of 

143 each analyzed field. The source detection criterion was based on a signal-to-noise ratio which is 

144 defined, assuming a Poisson distribution, as a ratio of the excess events above a background to its 

145 standard deviation. Photon counts were derived from each source and background regions and we 

146 set the detection threshold at 4(7. The source positions and the corresponding errors were obtained 

147 by fitting a 2D Gaussian around each X-ray spot. The source detection results are summarized in 

148 Tabled 



149 The light curves were constructed for each potential X-ray counterpart of the observed Fermi 

150 objects. Each light curve provides net-counting rates, with the count rates of the corresponding 

151 background region subtracted. In the timing analysis, the FI (XISO, XIS3) and BI (XIS1) CCD's 

152 light curves were combined using lcmath, and then re-binned using lcurve. To assess statisti- 

153 cal significances of the flux variations, the x 2 test was applied to each constructed dataset (probing 

154 a constant flux hypothesis with lest at s command). Finally, the XIS spectra for each source 

155 region were extracted, with the same corresponding background spectra as defined in the image 

156 analysis (see above). RMF files for the detector respo nse and ARF files fo r the effective area were 

157 generated using xisrmfgen and xissimarfgen dlshisaki et al.ll2007i) . In this spectral analy- 

158 sis, all the selected data from the FI CCDs were co-added (using mathpha) without calculating 

159 Poisson errors, and the response files were combined with the mart rmf and addrmf commands. 

160 Since all the studied FermifLAT objects are located at high Galactic latitudes, the absorption of 
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161 soft X-ray photons was se t to the Galactic one with the equivalent column density of a neutral 

162 hydrogen, iV H , as given in lDickey & Lockmanl (|1990|) . In some cases where apparent systematic 

163 features are visible as trends of the residuals with energy (see Figure [3]), we attempted to use an 

164 inter-calibration constant between the FI and BI CCDs to improve the fits. From this inspection, 

165 we found negligible improvement of the fits thus we conclude that the limited photon statistics is 

166 the predominantly responsible for the somewhat unsatisfactory model fits to the data. 



167 3. Results 

,68 3.1. 1FGLJ1231.4-1410 

169 Our Suzaku observations (interrupted for ~ 20 days3) revealed one X-ray point source (RA, 

170 Dec) = (187.°790(1), -14.°192(1)) within the LAT error circle of 1FGLJ123 1.4-1410. Figured] 

171 shows the corresponding X-ray image, prepared as described in § |2.2[ For further analysis, the 

172 source extraction region was set to within a 3' radius around the X-ray intensity peak, and the 

173 corresponding background region was chosen as indicated in Figure [TJ The light curve of the X- 

174 ray source with a time bin of 900 s is presented in Figure [2] The upper panel shows the count rate 

175 variation during the 1st observation, while the bottom panel shows that of the 2nd observation. The 

176 light curves of the two periods can both be well fitted by a constant count rate of 3.03 x 10~ 2 ct s _1 , 

177 with x 2 /d.o.f. = 58.3/107. This indicates that the X-ray emission of the analyzed source is steady, 

178 with the x 2 probability for a constant flux > 99%. 

179 The X-ray spectrum of the Suzaku source, which we propose to be the most likely counterpart 

180 of 1FGLJ1231. 1-1410, is shown in Figure [3j The energy range used for the fitting was set as 

181 0.4 — 7.0 keV. Fir st, we fit the X-ray spectrum by a blackbody emission moderated by the Galactic 

182 absorption only ( Morrison & McCammonlll983l) . This fit was not acceptable, however, due to 



183 significant residuals above 2keV (x 2 /d.o.f. = 128.1/34, see Figure[3l left panel, where the excess 

184 emission above 2keV has been enhanced by fixing the black body parameters to those determined 

185 by the data below 2 keV only). The situation was essentially unchanged when the column density 

186 was treated as a free parameter. In this case, the residuals above 2keV remained, but the returned 

187 value of Nh was then consistent with zero. To account for the > 2keV emission, we therefore 

188 added a power-law component to the thermal one, and fixed N a = 0. The goodness of the fit was 

189 in this way substantially improved to x 2 of 55.46/32, supporting the presence of a non-thermal tail 

190 in the spectrum of the analyzed object (see Figure[3]). In order to further confirm the reality of this 

191 finding, we analyzed the highest quality FI CCD (XIS0, XIS3) data which had sufficient photon 



'The exposure was interrupted because of the Target of Opportunity observation of GRB 090708. 
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192 statistics within the 2 — 8 keV range, examining various approaches for the background estimation, 

193 namely (i) the background taken from the same CCD chips but off-axis, as given in Figured! (ii) 

194 the concentric ring background surrounding the source region on the same CCD chips, and (iii) the 

195 background for the same region as the source estimated from the Lockman Hole observation taken 

196 with the same XIS mode at nearby dates (OBS ID = 104002010). In all of the examined approaches 

197 the presence of the non-thermal component in the X-ray spectrum of 1FGLJ1231. 1-1410 could 

198 be confirmed at high significance, as summarized in Tables |4] and [5J 

199 To sum up, we conclude that the X-ray counterpart of 1FGLJ1 23 1.1-1410 is robustly charac- 

200 terized by a blackbody-type spectrum with a temperature of kT ~ 0.16 ± 0.03 keV plus a power- 

201 law tail with the photon index of T ~ 1.79lo ^. The energy flux of the non-thermal component 

202 is 5*2-8 keV — 5.81 x 10~ 14 ergcm~ 2 s -1 , which can be compared with the FermifLAT energy flux 

203 So.i_ioGeV — (1-06 =t 0.06) x 10~ 10 erg cm -2 s~\ as given in the 1FGL catalog. Thus, the extrap- 

204 olation of the X-ray power-law component to the 7-ray range yielding the 0.1 — lOGeV energy 

205 flux ~ 5.74 x 10~ 13 erg cm~ 2 s -1 , falls orders of magnitudes below the observed one. This implies 

206 either a multi-component character or a concave spectral form of the high-energy X-ray-to-7-ray 

207 continuum of 1FGLJ123 1.1-1410, and both possibilities should be kept in mind in the context of 

208 a very likely association of the discussed source with a MSP . Indeed, the MSP PSR J1231-1411 

209 (marked by a white cross in FigureQ]) was recently found by lRansom et al. ( 2010 ) via the detection 

210 of radio pulsations with the pulse period of 3.68 ms within the LAT error circle of 1FGL J1231.1- 

211 1410 using the Green Bank Telescope (GBT), just after our Suzaku observations. In addition, the 

212 Fermi spectrum shows a cut-off at a round a few GeV, which is consistent with the typical spec- 

213 trum of MSPs ([Ransom et alj|2010|) . The X-ray emitter obs erved by Suzaku is located roughly 

214 40" away from the newly discovered MSP PSR J123 1-141 1 (IRansom etalJboiol see Figure [B, 

215 but considering the limited pointing accuracy of the Suzaku/XIS (< 1'), both objects can be con- 

216 sidered as co-spatial. In fact, as described in IRansom et al.l (|2010|) . a SwiftfXRH source at (RA, 

217 Dec) = (187.7972,-14.1953) coinciding with the Suzaku one, was found to be positionally con- 

218 sistent (within the 90% error of 5."5) with that of the MSP PSR J123 1-141 1. 



219 3.2. 1FGLJ1311.7-3429 

220 Two X-ray point sources were found within the LAT error circle of 1FGLJ131 1.7-3429. 

221 Figure 0] shows the corresponding X-ray image with the northern Suzaku object, src A, located 

222 at (RA, Dec) = (197.°939(l), -34.°508(2)) and the southern source, srcB, positioned at (RA, 

22 3 Dec) = (197.°942(l), -34.°534(2)). Note that srcB is situated just marginally within the edge 

224 of the FermifL AT error circle. For the further analysis, we set the source regions to within V radii 

225 around the respective X-ray flux maxima. The derived light curves of src A and src B with time 
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226 bins of 450 s are presented in Figure [5] (upper and lower panels, respectively). As shown, during 

227 the first 20 ksec of the observation, src A exhibited a very rapid X-ray flare, with the count rate 

228 changing by a factor of 10. After the flare, however, the X-ray flux of src A remained steady. A 

229 constant fit to the light curve of src A returns x 2 /d.o.f. = 403.9/97, and hence the variability of 

230 this source is statistically significant. On the other hand, src B was characterized by a constant flux 

231 over the duration of the exposure (% 2 /d.o.f. = 45.0/97) with a count rate of 1.3 x 10~ 2 ct s -1 . 

232 Figure [6] shows the spectra of src A and srcB within the energy range 0.4 — 8.0keV. The 

233 best model fits for both newly discovered X-ray objects consist of power-law continua with photon 

234 indices T ~ 1.38 ±0.13 (src A) and T ~ 1.34±0.16 (src B), moderated by the Galactic absorption. 

235 The detail of the model fitting are summarized in Table [6] Note that the observed X-ray spectra 

236 of the two sources are very similar, and the X-ray fluxes of the objects are almost identical. It 

237 is important to emphasize at this point that because of the relatively large PSF of Suzaku/XIS (a 

238 half power diameter of ~ 3'), it is quite difficult to separate completely src A and src B — located 

239 only 1/6 apart — for the purpose of the spectral analysis. As a result, even though it is clear we 

240 are dealing with two physically distinct X-ray sources (each detected at high significance), their 

241 spectral parameters cannot be accessed robustly. 



242 3.3. 1FGLJ1333.2+5056 

243 Our Suzaku observations revealed multiple regions of enhanced X-ray emission inside the 

244 LAT error circle of 1FGLJ1333. 2+5056, as shown in the corresponding X-ray image in Figure [7J 

245 The associations of these faint X-ray sources with 1FGL J1333. 2+5056 are therefore quite ambigu- 

246 ous. Within the FermifLAT error circle covered by the XIS exposure^, five X-ray enhancements 

247 have been found with detection significances of more than 4er, and these are denoted here as src A, 

248 B, C, D and E (see Figure[7Jand Table©. 

249 The light curves of src A, B, C, D and E with 5760 s binning are shown in Figure [8] in the 

250 descending order. As noted above, all the analyzed X-ray sources are very dim, with X-ray fluxes 

251 at the level of ~ 10~ 14 ergcm~ 2 s -1 . Hence, we could not assess the variability properties of the 

252 selected objects by means of the \ 2 test with a constant flux hypothesis (see Table UJ. The spectra 

253 of the five X-ray sources, all extracted within 1' source radii, are shown in Figured Again, limited 

254 photon statistics precluded any detailed analysis, and therefore in the model fitting we applied 

255 only single power-law models moderated by the Galactic absorption. The results are summarized 



2 Note that the 1FGL localization error for the analyzed 7-ray object did not improve sufficiently between 0FGL 
and 1FGL. For this reason, we could not cover the entire 95% LAT error circle of 1FGLJ1333. 2+5056 within one 
pointing of Suzaku/XIS. 
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256 in Table [8] We also emphasize that the 1FGL error circle unfortunately runs off the edge of Suzaku 

257 field of view. For all these reasons, we cannot persuasively identify an X-ray counterpart of the 7- 

258 ray source 1FGLJ1333. 2+5056. Nevertheless, we note that one of the X-ray enhancements, src D, 

259 c oincides with th e z = 1.362 FSRQ CLASS J 1333+5057 (marked in Figure [7] by a white cross; 



Shaw et all 120091) . listed in the 1FGL as a possible association with 1FGLJ1333. 2+5056. Note 



261 however a relatively low significance of the detection of this source with Suzaku/XIS. 



3.4. 1FGLJ2017.3+0603 



263 A single prominent X-ray point source was found at the edge of the 1FGL error circle of 

264 the unidentified 7-ray source 1FGLJ2017. 3+0603. The X-ray source is located at (RA, Dec) = 

265 (304. °3 10(1), 6.°052(1)), as shown in Figure [TOl For the further analysis, we set the extraction 

266 region to encircle this bright source with a radius of 3'. The corresponding light curve of the newly 

267 discovered X-ray source is show in Figure QT] with 620 s binning. The light curve is consistent (at 

268 the level of > 99%) with a constant X-ray flux within the errors (x 2 /d.o.f. = 26.4/56) and the 

269 average count rate 4.07x 10~ 2 ct s -1 . Figured!] shows the X-ray spectrum of the analyzed source. 

270 A power-law model (photon index T ~ 1.6) with the Galactic absorption provided the best fit to 

271 the data, and the obtained best fit parameters are given in Table [9] 

272 The X-ray point source found at the edge of the 1FGL error circle is positionally coincident 

273 (offset by 15", which is muc h less the Suzaku /XIS position accuracy of ~ 1') with the z = 1.743 

274 FSRQ CLASS J2017+0603 dMvers et al.l l2003h. This blazar has been listed in the first Fermi/LAI 

275 AGN Catalog (Ubdo et alJboiOdh as being possibly associated with 1 FGL J20 1 7 .3+0603 , even 

276 though the probability for such an association was not quantified. We denote its position in Fig- 

277 ure[10]with a white cross. More recently, radio and 7-ray pulsations with the pulse period of 2.9 ms 

278 have been discovered using the Nancay rad io telescope well within the Fermi/LAtI error circle of 

279 1FGL J2017. 3+0603 (|Cognard et al.ll2010|) . pointing instead to a pulsar (rather than blazar) as- 

280 sociation of this source. In Figure [10] we mark the position of the MSP PSR J2017+0603 with 

281 a black cross. As shown, no X-ray counterpart of the pulsar has been detected by Suzaku/XIS. 

282 In order to determine the corresponding X-ray flux upper limit, we set an additional source re- 

283 gion within 1' radius around the position of the radio pulsar, and assumed a power-law emis- 

284 sion spectrum with photon index T = 2. The resulting 90% confidence X-ray upper limit is 



S- 



2-8 kcV 



< 2.61 x 10- 14 ergcm- 2 s- ] 
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286 



4. Discussion 



4.1. The Observed Sample 



288 Within the error circle of the unidentified 7-ray object 1FGL J1231.4— 1410, only one X-ray 

289 source was found, w hich is positionally consistent with the radio/7-ray MSP PSR J123 1—141 1 

290 dRansom et alJboid see Figured]). The broadband spectrum of 1FGL J1231.1-1410/PSRJ1231- 

291 1411, including our Suzaku/XIS data and the derived UVOT optical/UV upper limits from Swift, 

292 are shown in Figure [131 We note that the SE P is reminiscent of that of the Geminga pulsar 

293 ([Thompson et al.lll999l) . or 3EG J1835+5918 (IHalpern et all 120021) . Hence the consistent picture 

294 emerges, in which the kT ~ 0. 16 keV blackbody component of the newly discovered X-ray coun- 

295 terpart of 1FGLJ123 1.1-1410 originates as thermal emission from the surface of a rotating mag- 

296 netized neutron star, a pulsar, while the 7-ray emission detected by Fermi/LAtI may be accounted 

297 by the emission of ultra-relativistic electrons present within the pulsar magneto sphere. The non- 
298 thermal X-ray component is then likely to be produced within the magnetosphere of PSR J 1231- 
299 1411 as well, even though one may also expect some contribution from the outer regions (pulsar 

300 wind nebulae) to the detected hard X-ray continuum. 

301 Assuming that PSR J1231— 141 1 is a typical MSP with th e spin period P = 3.68 ms and a 

302 spin-down rate P = 2.1 X 10~ 20 s s _1 (see lRansom et al.ll2010h . one can calculate the correspond- 

303 ing spin-down luminosity as L S( j = A-k 2 IPP^ ~ 2 x 10 34 ergs _1 , and the magnetic field in- 

304 tensity at the light cylinder (radius, R = cP/2n) as B ic ~ 4vr 2 (3LP/2c 3 P 5 ) 1/2 ~ 5 x 10 4 G, 

305 where / = 1 45 g cm 2 is the sta r's moment of inertia. Meanwhile, for the claimed distance 

306 d ~ 0.4kpc (|Ransom et al.ll2010|) . the observed 7-ray luminosity of PSR J1231— 141 1 leads as 



307 L 7 ~ 2 x 10 33 ergs~ 



its non-thermal X-ray luminosity is Lx — 10 ergs , and the total X- 



308 ray luminosity Lx/tot ~ 3 x 10 ergs . These values are then consistent with the millisec- 

309 ond pulsar scenario - outer- magnetosphere mod els in particular - in a framework of which one 

310 should expect L 7 ~ 0.1 L.a (lAbdo et alJl2009dh and L x ~ 10~ 3 L sd feecker & Truemper![l997 : 



311 



Gaensler & Slanel I2006I : IZhang et al.l |2007|) . with relatively large dispersion, however. Interest- 



312 ingly, the synchrotron X-ray luminosity produced close to the light cylinder within the expected 

313 magnetic field B\ c and a fraction (say, 10%) of the volume V ~ -R 3 , would then be close to 

314 the observed non-thermal X-ray luminosity assuming rough energy equipartition between ultra- 

315 relativistic electrons and the magnetic field. 

316 In the case of 1FGLJ1 31 1.7-3429, two potential X-ray counterparts have been discovered 

317 in our Suzaku observations. The association of this Fermi object with the northern source srcA 

318 is more likely, since the southern X-ray spot srcB is located only marginally within the 95% 

319 FermifL AT error circle of the 7-ray emitter. Yet the classification of 1FGL J131 1.7-3429/src A, 

320 for which the broad-band spectrum (including radio and optical upper limits) is shown in Fig- 
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321 ure[[4l remains vague. Currently, no radio or 7-ray pulsations have been found at the position of 

322 1FGLJ13 11.7-3429, and this favors an extragalactic origin of the detected high-energy emission. 

323 And indeed, the flat X-ray continuum (r ~ 1.4) and the 7-ray-to-X-ray energy flux ratio > 100 

324 (with So.i_iogcV — 6.4 x 10~ n ergcm~ 2 s _1 as given in the 1FGL catalog) would be c onsistent 

325 with the characteristics of luminous blazars of the FSRQ type (e.g., ISikora et al.ll2009|) . On the 

326 other hand, however, the radio upper limit indicating the GHz energy flux ~ 10~ 5 times smaller 

327 than the GeV energy flux, invalidates the blazar nature of 1FGLJ1311. 7-3429. That is because 

328 all active galaxies established till now as 7-ray emitters are characterized by relatively strong, 

329 Doppler-boosted radio emission. In particular, radio energy fluxes of bona fide blazars included 

330 in OFGL are, for a given FermifLAT photon flux of ~ 10~ 7 photons cm~ 2 s _1 , at least an order 

331 of magnitude highe r than the energy flux implied by the NVSS upper limits for src A (see, e.g., 



332 



Kovalev et al.ll2009|) . In addition, a very prominent lOks-long X-ray flare detected from src A, to- 



333 gether with the steady GeV flux of 1 FGL J 1 3 1 1 .7-3429, would not match easily a typical behavior 

334 of FSRQs: this class of blazars is known for displaying dramatic variability at GeV photon ener- 

335 gies, but only modest variations in the X-ray band. Therefore, the nature of the analyzed Fermi 

336 source and its newly discovered Suzaku counterpart remains an open question. 

337 Within the error circle of 1FGL J 1333. 2+5056, our Suzaku/XIS observations revealed the 

338 presence of several weak X-ray flux maxima with possibly diverse spectral properties (as indi- 

339 cated by the spectral analysis hampered by the limited photon statistics). One of the detected 

340 X-ray sources (srcD) coincides with the high-redshift blazar CLASS J1333+5056 (z = 1.362). 
34, The broad-band spectral energy distribution of 1FGL J1333.2+5056/CLASS J1333+5056/srcD is 

342 presented in Figure [151 including the LAT 7-ray, Suzaku X-ray, archival radio, and newly ana- 

343 lyzed Swift/UVOT data for the blazar. The constructed SED reveals two distinct radiative compo- 

344 nents, consisting of a low-energy synchrotron bump and an (energetically dominant) high-energy 

345 inver se-Compton continuum , reminiscent of typical broad-band spectra for blazars of the FSRQ 

346 type ( Ghisellini et al. 1998|) . Note that the X-ray-to-7-ray flux ratio ~ 10 3 implied by Figure [T5l 

347 as well as the relatively large radio flux, would be both in agreement with the blazar identification 

348 of 1FGLJ1333. 2+5056. In addition, we note that the discussed Fermi object is the most variable 

349 in 7-ray s out of all four Fermi targets studied in th is paper, with the v ariability index of 38 (which 

350 indicates a < 1% probability of a steady flux; see lAbdo et al.ll2010a) . The additional support for 

351 the blazar association is offered by the fact that the 7-ray continuum of 1FGLJ1333. 2+5056 is the 

352 steepest among the four Fermi objects observed by us, with the photon index ~ 2.5 ± 0.1, which 

353 is compatible with t he mean 7-ray ph oton index of the FSRQ population reported in the 1FGL, 

354 namely 2.47 ± 0.19 JAbdo et allboipj) . 



355 Finally, in the case of 1FGL J20 17.3+0603, the MSP PSR J2017+0603 was newly discov- 

356 ered by the Nancay Radio Telescope (ICognard et alj 120101) . and the association between the ra- 

357 dio and 7-ray sources was confirmed by the pulse detection with the same period in the LAT 
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358 data. Interestingly, in our Suzaku/XIS exposure we have only detected the high-redshift blazar 

359 (z = 1.74 3) CLASS J20 17 +0603, but not the pulsar. The same is true for the SwiftflTVGT ob- 

360 servation (|Cognard et al.ll2010|) . which resulted in analogous flux and upper limit measurements 

361 in the optical for the blazar and pulsar, respectively. The constructed radio to X-ray SEDs for the 

362 pulsar and blazar systems ar e shown in Figure [16] together with the LAT spectrum. Regarding the 

363 pulsar. ICognard et alJ(|2010|) discovered that PSR J2017+0603 is located at a distance d ~ 1.6 kpc, 

364 and as such is characterized by the spin-down luminosity L sd ~ 1.34 x 10 34 ergs _1 . The X-ray 

365 (2 — 8keV) luminosity derived from the Suzaku/XIS upper limit for this pulsar, L x < 8.0 x 

366 1 30 erg s _1 , is then consistent with the expected "pulsar-like" luminosity ratio Lx/^sd ~ 10~ 3 . 

367 The overall curved 7-ray spectrum of 1FGLJ2017. 3+0603, characterized by the small photon 

368 index ~ 1.88 ± 0.05, supports the pulsar association. On the other hand, the relatively large 

369 radio flux of CLASS J20 17 +0603, together with the X-ray-to-7-ray flux ratio ~ 300 for the 

370 1FGLJ2017.3+0603/CLASS J2017+0603 system, are in some level of agreement with the blazar 

371 interpretation. The 7-ray photon index of 1FGLJ2017. 3+0603 is however rather flat for a FSRQ 

372 and represents a ~ 3cr deviation from the distribution observed for FSRQs (mean= 2.47, a = 0.19; 



373 see 



Abdo et al.ll2010fh thus making the association with the FSRQ less likely. Although the de- 



374 tected pulsations in radio and 7-rays is key to the identification of the 7-ray source with a pulsar, 

375 there may be some contaminating flux from the blazar. Indeed, the chance probability of finding a 

376 CLASS-like background blazar in the Fermi error circle of this source is ~ 0.003%. Considering 

377 over 1400 sources in the 1FGL catalog, such 'mixed' cases could be expected. 



378 



4.2. Implications 



379 What class of astrophysical objects can be in general associated with the unidentified high 

380 Galactic latitude 7-ray sources? It was noted, for example, that compact and relat ively nearby 



381 molec ular clouds exist at \b\ > 10°, and these should emit 7-rays at least at some level. iTorres et al 



382 (|2005|) argued, however, that the expected GeV emission of such clouds is too low to account for 

383 the observed fluxes of unidentified EGRET sources, and the same applies to the bright unidentified 

384 Fermi/LAT objects. Anoth er classes of possible counterpa rts proposed were radio-quiet pulsars and 

385 isolated neutron stars (e.g. jYadigaroglu & Romarii 1995 ). and this idea has indeed been validated 

386 by the subsequent multi-frequency studies, as discussion in § [Q We note in this context that the 

387 Galactic origin of high-latitude 7-ray emitters is especially probable for the objects located at 10° < 

388 \b\ < 30° within the Gould Belt (~ 0.3 kpc from the Earth), whic h constitutes an agg regation of 

389 massive late-type stars, molecular clouds, and supernova remnants (|Grenier et al. 2000). 



390 A probably more challenging population of 7-ray emitters is represented by the isotropic com- 

391 ponent of the unidentified EGRET objects, consisting of about 60 sources (about one third of which 
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392 with the G alactic latitudes \b\ > 45°, including several n on/weakly-variable duri ng the EGRET ob 



393 servations: 



Ozel & Thompsonlll996l: Geh rels et al.l 120 00). For those sources. iTotani & Kitav ama 



394 (|2000|) have for example suggested associations with large-scale shocks prod uced during the struc- 

395 ture formation in the intergalactic medium (see also Iwaxman & Loeb 2000 ). Totani & Kitayama 

396 explored the connection between steady GeV objects located off the Galactic plane, and labeled 

397 in the 3EG catalog as "possibly extende d," with dynamically form ing clusters of galaxies (and 

398 not single virialized cluster systems; see iKawasaki & Totanill2002|) . However, the non-variable 



399 natur e of the 7-ray emission of several of the considered objects was questioned (see lReimer et al. 



400 



2003L and references therein), and the high efficiency of the particle acceleratio n at the structure 



401 formation shocks required by the model was also noted (e.g jKeshet et al.ll2003|) 



402 Radio galaxies are prime candidates for the unidentified high Galactic latitude EGRET sources, 

403 especially since the only confirmed non- blazar AGN detected previously at GeV photon energies 

404 was the nearby radio galaxy Centaurus A (lSteinle et al.ll 19981 : ISreekumar et al.lll999l) . Yet no other 

405 radio galaxy has been firmly detected b y EGRET at the sign ificance leve l high enough (> 4cr) to 

406 be included in the 3rd EGRET catalog dHartman etaDl 1999b . Moreover. Icillis et al.1 hoo4) . who 

407 applied a stacking analysis of the EGRET data for a sample of the brightest and/or the closest radio 

408 and Seyfert galaxies, showed that 'no detection significance greater than 2 a has been found for any 

409 subclass, sorting parameter, or number of objects co-added.' Nevertheless. iMukherjee et al.l(l2002|) 

410 argued that the most likely counterpart to the unidentified EGRET source 3EG J1621+8203 is the 

411 bright radio galaxy NGC 6251. A m arginal detection o f 3C 1 1 1 with EGRE T has also been reported 

412 (|Hartman. Kadler & Tuellerll2008|) . We also note that ICombi et al.l (I2003 ) reported the discovery 

413 of a new radio galaxy within the location error circle of the unidentified 7-ray source 3EG J1735- 

414 1500. The identification of 3EG J1735-1500 was however controversial, due to the presence of an 



415 anoth er likely (blazar-type) candidate within the EGRET error contours (ISowards-Emmerd et al. 



2004k . The most recent analysis based on the 1 5 months of FermilhPS data resulted in the detection 



417 of 1 1 non-b lazar-type AGN (a ll radio galaxies), including the aforementioned cases of NGC 6251 

418 and 3C1 1 1 (|Abdo et al.ll2010gh . The idea that some fraction of unidentified 7-ray emitters may be 

419 associated with faint radio galaxies is therefore validated, although this should rather apply to only 

420 dimmer Fermi objects, and not to the population of exceptionally bright 7-ray sources detected 

421 already by EGRET. 

422 The Suzaku/XIS studies of four bright FermifLAT objects reported here provide an impor- 

423 tant contribution to the debate regarding the nature of unidentified 7-ray emitters located at high 

424 Galactic latitudes. In particular, our observations support the idea that a significant fraction of 

425 such objects may be associated with old (> Gyr) MSPs present within the Galactic halo and the 

426 Earth's neighborhood (such as 1FGLJ1231. 1-1410 and 1FGLJ2017. 3+0603). Yet not all of the 

427 unidentified Fermi objects are related to the pulsar phenomenon. Instead, some of those may 

428 be hosted by active galaxies, most likely by the luminous and high-redshift blazars of the FSRQ 
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429 type (1FGLJ1333. 2+5056 is as good blazar candidate, for example). However, there still remain 

430 unidentified sources, (e.g., 1FGLJ131 1.7-3429), for which neither blazar nor pulsar scenarios 

431 seem to apply. For these, ultra-deep multi-wavelength studies are probably needed to unravel their 

432 physical nature. 



433 5. Summary 

434 In this paper we reported on the results of deep X-ray follow-up observations of four uniden- 

435 tified 7-ray sources detected by the Fermi fLAT instrument at high Galactic latitudes (\b\ > 10°) 

436 using the X-ray Imaging Spectrometers onboard Suzaku satellite. All of the studied objects have 

437 been detected at high significance (> lOer) during the first 3-months of the Fermz'/LAT operation. 

438 For some of them, possible associations with pulsars and blazar have been recently discussed, and 

439 our observations provide an important contribution to this debate. In particular, an X-ray point 

440 source was newly found within 95% error circle of 1FGLJ1231. 1-1410. The X-ray spectrum of 

441 the discovered Suzaku counterpart of 1FGLJ1231. 1-1410 is well fitted by a blackbody emission 

442 with a temperature of kT ~ 0.16keV plus an additional power-law component with a differential 

443 photon index T ~ 1.8. This supports the recently claimed identification of this source with a MSP 

444 PSRJ123 1-141 1. For the remaining three Fermi objects, the performed X-ray observations are 

445 less conclusive. In the case of 1FGLJ1311. 7-3429, two possibly associated X-ray point sources 

446 were newly found. Even though the 0.4 — lOkeV spectral and variability properties for those 

447 could be robustly accessed, the physical nature of the X-ray emitters and their relations with the 

448 7-ray source remain unidentified. Similarly, we found several weak X-ray sources in the field of 

449 1FGLJ1333.2+5056, one coinciding with the high-redshift blazar CLASS J1333+5057. We ar- 

450 gue that the available data are consistent with the physical association between these two objects, 

451 even though we were not able to identify robustly the Suzaku counterpart of 7-ray emitter due to a 

452 large positional uncertainty of 1FGLJ1333. 2+5056. Finally, we found an X-ray point source in the 

453 vicinity of 1FGLJ2017. 3+0603. This Fermi object was recently suggested to be associated with 

454 a newly discovered MSP PSR J2017+0603 because of the detection of radio and 7-ray pulsations. 

455 However, we did not detect the X-ray counterpart of the pulsar, but instead of the high-redshift 

456 blazar CLASS J2017+0603 located within the error circle of 1 FGL J20 1 7 .3+0603 . Still, the result- 

457 ing upper limits for the X-ray emission do not invalidate the pulsar association. 



458 
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Table 1 : EGRET and FermifLAI entries for the analyzed objects 



Name 


RA 


DEC 


I 


b 


-fo.l-20 GeV 


^95% 




[deg] 


[deg] 


[deg] 


[deg] 


[10 _8 phcm- 2 s _1 ] 


[deg] 


1FGLJ1231. 1-1410* 


187.80 


-14.17 


295.53 


+48.41 


14.9+0.7 


0.03 


3EGJ1234-1318 


188.19 


-16.30 


296.43 


+49.34 


21.6+5.3 


0.76 


(EGRJ123 1-1412) 














lFGLJ1311.7-3429§ 


197.95 


-34.49 


307.69 


+28.19 


11.7+1.1 


0.04 


3EGJ1314-3431 


198.51 


-34.52 


308.21 


+28.12 


18.7+3.1 


0.56 


(EGRJ1314-3417) 














lFGLJ1333.2+5056§ 


203.30 


+50.94 


107.32 


+64.90 


4.5+1.0 


0.15 


3EGJ1 337+5029 


204.39 


+50.49 


105.40 


+65.04 


9.2+2.6 


0.72 


(EGRJ1338+5102) 














1FGLJ2017.3+0603* 


304.34 


+6.05 


48.62 


-16.02 


4.5+0.5 


0.04 



' Data consistent with no variability between EGRET and Fermi/LAI observations. 

§ The 7-ray fluxes determined by EGRET and Fermi/LAI marginally consistent within 2a level. 

■t Data consistent with no variability between EGRET and Fermi/LAI observations because of the EGRET 

detection limit ~ 6 x 10~ 8 phcm -2 s _1 . 
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Table 2: Suzaku/XIS Observation Log 



Name 


OBS ID 


Pointing 


Center* 


Observation start 


Effective exposure 






RA [deg] 


DEC [deg] 


(UT) 


[ksec] 


1FGLJ1231. 1-1410 


804017010 t 


187.8001 


-14.1665 


2009/07/08 22:53:48 


23.8 




804017020 t 


187.7993 


-14.1672 


2009/07/28 05:21:37 


44.8 


1FGLJ1311.7-3429 


804018010 


197.9603 


-34.4918 


2009/08/04 04:56:35 


33.0 


1FGLJ1333.2+5056 


804019010 


203.2955 


51.0170 


2009/06/01 10:13:15 


39.1 


1FGLJ2017.3+0603 


804020010 


304.3461 


6.0496 


2009/10/27 10:14:45 


26.7 



* The pointing centers were the positions given in 0FGL catalog (lAbdo et al.ll2009ch . 
t The requested continuous 80 ks Suzaku exposure was interrupted by Target of Opportunity (ToO) observa- 
tion of GRB 090708. For this reason the observation was divided into 30 ks and 50 ks segments before and 
after the ToO observation. 



Table 3: Source detection results of Suzaku observation 



Name Position Detection Significance r 95 % 







RA [deg] 


DEC [deg] 


a 


[arcsec] 


1FGLJ1231.4-1410 




187.790 


-14.192 


13.67 


7.44 


1FGLJ1311.7-3429 


src A 


197.939 


-34.508 


15.52 


17.44 




srcB 


197.942 


-34.534 


12.89 


12.27 


1FGLJ1333.2+5056 


src A 


203.252 


50.983 


8.53 


23.34 




srcB 


203.161 


51.032 


7.27 


19.97 




srcC 


203.276 


51.014 


7.47 


20.75 




src D 


203.479 


50.967 


4.50 


38.41 




srcE 


203.381 


50.892 


4.91 


26.51 


1FGLJ2017.3+0603 




304.310 


6.052 


14.44 


4.73 
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Table4: Fitting Parameters for 1FGL J 123 1.1 — 1410 in a framework of blackbody (BB) and power- 
law (PL) models 





BB model 


BB+PL Model 


parameter 


value & error 


value & error 


N H [10 22 cm" 2 ] 


0.0 (fixed) 


0.0 (fixed) 


kT [keV] 


0.228 ± 0.008 


lfi o +0.024 
U.lOO _ .026 


norm. (BB) 


(1.42 ±0.14) x 10~ 6 


(1.20 t° f 7 ) x 10-6 


r 




1 70 +0-40 
L - ' y -0.17 


norm. (PL) 




(1.94 x 10- 5 


x 2 


128.1 


55.46 


d.o.f. 


34 


32 


reduced \ 2 


3.768 


1.733 


Flux (2 -8keV) 




(5.79 ±{f 2 ) xlO- 14 



[erg cm 2 s 1 ] 
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Table 5: Blackbody (BB) and power-law (PL) components in the X-ray spectrum of 
1FGLJ1231. 1-1410 





(i) Standard Background 


(ii) Ring Background 


(iii) Lockman Hole Background 




BB 


BB+PL 


BB BB+PL 


BB 


BB+PL 


x 2 


134.71 


56.16 


67.21 18.97 


39.77 


23.97 


d.o.f. 


34 


32 


30 28 


38 


36 


F value 




22.4 


35.6 




11.9 


Probability 


8 


.33xl0" 7 % 


2.04 xlO" 6 % 




LlOxlO" 2 % 
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Table 6: Fitting Parameters for 1FGL J1311.7— 3429 for power-law model 





src A 


srcB 


parameter 


value & error 


value & error 


N H [10 20 cm" 2 ] 

r 

norm. 


4.45 (fixed) 

i oo +0.13 
L - 6ii -0.13 

(2.69 t° f 7 ) x 10~ 5 


4.45 (fixed) 

1 U +0.16 
L - oqi -0.15 

(2.08 toil) x 10" 5 


x 2 


42.6 


42.1 


d.o.f. 


38 


38 


reduced \ 2 


1.12 


1.11 


Flux (2-8keV) 
[ ergcm~ 2 s _1 ] 


(1.45 t° il) x lO" 13 


(1.20 x lO" 13 



-25- 



Table 7: Count rates and constant flux fits for X-ray sources within the error circle of 
1FGLJ1333.2+5056 



Source 


Average count rate & Error 


X 2 /d.o.f. 


Prob. 




[ 10~ 3 cts" 1 ] 




[%] 


src A 


5.47 ±0.51 


14.9/15 


46.08 


srcB 


4.40 ± 0.49 


22.4/15 


9.73 


srcC 


4.37 ± 0.48 


18.8/15 


22.28 


srcD 


2.19 ±0.44 


20.5/15 


15.25 


srcE 


1.70 ± 0.44 


23.3/15 


7.86 
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Table 8: Fitting Parameters for 1FGLJ1333. 2+5056 for power-law model 





src A 


srcB 


srcC 


srcD 


srcE 


parameter 


value & error 


value & error 


value & error 


value & error 


value & error 


N H [10 20 cm- 2 ] 

r 

norm. [xlCT 5 ] 


1.09 (fixed) 

9 or +0.35 
-0.32 

1 ry+0.28 

l - 0< -0.28 


1.09 (fixed) 
1.48 i°; 29 7 

lot 18:25 


1.09 (fixed) 

1 r.1 +0.31 

L - 0L -0.29 

0.84 t° H 


1.09 (fixed) 
1 76 +0 - 60 

L - 10 -0.52 

77 +0 - 31 

u -' ' -0.30 


1.09 (fixed) 

9 r. n +0.69 
z - ou -0.58 

i +0 - 36 
J-" 54 -0.37 


x 2 


13.0 


7.33 


18.3 


12.7 


11.4 


d.o.f. 


18 


18 


18 


16 


12 


reduced \ 2 


0.720 


0.407 


1.02 


0.796 


0.949 


Flux (2-8keV) 
[xlO' 14 ergon^s -1 ] 


9 1fi +0.88 
z - iD -0.75 


4.98 ±{f 7 


o 77 +1.17 
' ' -1.11 


9 41 +1.55 
z - 4i -1.26 


-L-^ -0.94 
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Table 9: Fitting Parameters for 1FGLJ2017. 3+0603 for power-law model 



parameter value & error 

N H [10 22 cm- 2 ] 

r 

norm. 

X 2 34.8 
d.o.f. 38 
reduced % 2 0.916 

Flux (2 - 8keV) (1.99 x 10 _ 

[ ergcm~ 2 s^ 1 ] 



0.1309 (fixed) 



1.59 



+0.15 



(5.03 



0.15 

+0.68^ 1f) -5 
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Fig. 1.— Suzaku/XIS FI (XIS0+3) image of 1FGLJ1231. 1-1410 region in the 0.4 - lOkeV 
photon energy range. The image shows the relative excess of smoothed photon counts (arbitrary 
units indicated in the bottom bar) and is displayed with linear scaling. The areas enclosed by 
solid and dashed circles are source and background regions, respectively. Thick solid circle de- 
notes 95% positio n error of 1FGLJ123 1.1 — 1410. White cross marks the position of radio MSP 
PSRJ1231-1411 kansom et alhoioh . 
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Fig. 2. — Suzaku/XIS light curves of the X-ray counterpart of 1FGL J 123 1.1 — 1410 during the 1st 
and the 2nd observations (upper and lower panels, respectively). Binning time applied is 900 s. The 
zero point of the upper and lower panels are MJD 55020.9971 and 55040.2343 (TDB: Barycentric 
Dynamical Time). 
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(a) Blackbody Model (b) Blackbody+Power-Law Model 




Energy (keV) Energy (keV) 



Fig. 3. — Suzaku/XIS spectra of the X-ray counterpart of 1FGLJ1231. 1 — 1410 in the photon en- 
ergy range 0.4 — 7.0 keV fitted with the blackbody model (a) and blackbody+power-law model 
(b). FI data are shown in black, and BI data in gray. 
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Fig. 4.— Suzaku/XIS FI (XISO+3) image of 1FGL J13 1 1 .7-3429 region in the 0.4- 10 keV photon 
energy range. The image shows the relative excess of smoothed photon counts (arbitrary units 
indicated in the bottom bar) and is displayed with linear scaling. The regions enclosed by solid 
and dashed circles are source and background regions, respectively. Thick solid circle denotes 95% 
position error of 1FGLJ131 1.7— 3429. Within this error circle, two potential X-ray counterparts 
of the 7-ray source were found: src A and src B. 
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Bin Time: 450.0 s 




4x10 4 6x10 4 
Observation Time (s) 



Fig. 5. — Suzaku/XIS light curves of two potential X-ray counterparts of 1FGL J13 1 1 .7—3429 with 
450 s binning. The northern source src A (upper panel) showed highly significant X-ray flare in the 
first 20 ks of observation, during which the count rate increased by a factor of 10. The southern 
source src B (lower panel) was steady during the whole exposure. The zero point of src A and src B 
is MJD 55047.2280 (TDB). 
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Fig. 6. — Suzaku/XIS Spectra of two possible X-ray counterparts of 1FGLJ131 1.7— 3429 in the 
photon energy range 0.4 — 8.0 keV fitted with the best fit power-law model. FI data are represented 
in black, and BI data in gray. 
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i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




2E-08 4E-08 6E-08 8E-08 1E-07 1.2E-07 1.4E-07 

Fig. 7.— Suzaku/XIS FI (XISO+3) image of the 1FGLJ1333. 2+5056 region in the 0.4 - lOkeV 
photon energy range. The image shows the relative excess of smoothed photon counts (arbitrary 
units indicated in the bottom bar) and is displayed with linear scaling. The regions enclosed by 
solid and dashed circles are source and background regions, respectively. Thick solid ellipsis 
denotes 95% position error of 1FGLJ1333. 2+5056. Within this error circle, several potential X- 
ray counterparts of the 7-ray object were found. White cross marks the position of active galaxy 
CLASS J1333+5057. 
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Bin Time: 5760.0 s 
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Fig. 8. — SuzakufXIS light curves of potential X-ray counterparts of 1FGL J1333. 2+5056 with the 
applied time binning of 5760 s. The zero point of time is MJD 54983.4274 (TDB). 
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Fig. 9. — Suzaku/XIS spectra of the selected possible X-ray counterparts of 1FGLJ1333. 2+5056 
fitted with a power-law model. FI data are represented in black, and BI data in gray. 
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Fig. 10.— Suzaku/XIS FI (XIS0+3) image of the 1 FGL J20 17.3 +0603 region in the 0.4 - lOkeV 
photon energy range. The image shows the relative excess of smoothed photon counts (arbitrary 
units indicated in the bottom bar) and is displayed with linear scaling. The regions enclosed by 
solid and dashed circles are source and background regions, respectively. Thick solid circle denotes 
95% position error of 1FGLJ2017. 3+0603. One X-ray point source was found within this error 
circle. White cross mark denotes the position of the blaza r CLASS J2017+060 3. Black cross mark 
denotes the position of the radio MSP PSR J2017+0603 Jcognard et alJboioh . 
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Bin Time: 620.0 s 
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Fig. 11. — Suzaku/XIS light curve of an X-ray point source within the error circle of 
1FGLJ2017.3+0603 with the applied 620 s time binning. The zero point of time is MJD 
55131.4285 (TDB). 
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Fig. 12. — SuzakufXIS spectrum of the potential X-ray counterpart of 1FGLJ2017. 3+0603 with 
the best fit power-law model. FI data are shown black, and BI data in gray. 
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1FGLJ1231. 1-1410 




Energy [eV] 



Fig. 13.— Broad-band spectrum of 1FGL Jl 23 1.1- 1410/PSRJ1 23 1-1411. The X-ray data 
points represent the weighted mean of Suzaku/X IS FI and BI data. The 7-ray data points are 
taken from the 1FGL catalog (|Abdo et alJl2010q) . The rad io data point i s deriv ed from the MSP 
PSR Jl 23 1-141 1 observed with Green Bank Telescope bv lRansom et alJ ( bold) . The optical/UV 
upper limits were derived from the SwiftfUVOT observation 
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Fig. 14. — Broad-band spectrum of 1FGLJ131 1.7— 3429. The X-ray data points represent the 
weighted mean of Suzaku/XIS FI an d BI data for srcA. The 7-ray data points are taken from 
the 1FGL catalog ( Abdo et al.ll2010cr) . The radio upper limit is taken from the NVSS catalog 
JcondonetalJll998h . The optical/UV data points show the SwiftfUYOT data. 
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Fig. 15.— Broad-band spectrum of 1FGL J1333.2+5056/CLASS J1333+5057. The X-ray data 
points represent the weighted mean of Suzaku/XIS FI and BI data for src D w hich coincides wit h 
the CLASS source. The 7-ray data points are taken from the 1FGL catalog (lAbdo et al.ll2010ch . 
The radio dat a points, represent ing blazar CLASS J1333+505 7, are taken from the CLASS catalog 
(filled circle: iMvers et alJl2003h. N VSS catalog (filled square: bondon et alJll998h and GB6 cata- 
log (filled triangle; Gregory et al. 19961) . Optic al data point (open circle) was de rived from SDSS 
J133353.78+505735.9 (SDSS Data Release 6: kdelman-McCarthv et aDl2008h . optical/UV data 
points and upper limit (filled circle) is the Swift/UVOT data. 
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Fig. 16. — Broad-band spectrum of 1FGLJ2017. 3+0603. The X-ray data points represent the 
weighted mean of Suzaku/XIS FI and BI data for active galaxy CLASS J2017+0603. The X-ray 
upper limit (open square) is derived from t he location of the M SP PSR J2017+0603. The 7-ray data 
points are taken from the 1FGL catalog (|Abdo et al.l 1201 Pel) . The radi o data points, repr esenting 
CLASS J2017+0603, are taken from the C LASS catalog (filled circle: iMvers et al.ll2003h • NVSS 
catalog (filled square; ICondon et al] [T998) and GB6 catalog (filled triangle; Gregory et al.l ll996). 
The open diamon d shaped point in ra dio shows the MSP PSR J2017+0603 observed with Nancay 
Radio Telescope (ICognard et alJl2010l) and also the op tical/UV upper limits (open circle) show the 
MSP observed with Swift/UYOT ( Cognard et al.ll2010 ). The optical/UV data points (corresponding 
filled circle) show the blazar CLASS J2017+0603 observed with Swift/XJVOT. 



